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Abstract 

We present some applications of general harmonic/wavelet 
analysis approach (generalized coherent states, wavelet pa- 
ckets) to numerical/analytical calculations in (nonlinear) 
quasiclassical/quantum beam dynamics problems. (Naive) 
deformation quantization, multiresolution representations 
and Wigner transform are the key points. 



resembles the action of a microscope. We have contribu- 
tion to final result from each scale of resolution from the 
whole infinite scale of spaces. Besides affine group sym- 
metry, in part 3 we consider modelling, based on very use- 
ful and quantum oriented Wigner transform/function ap- 
proach (corresponding to Weyl-Heisenberg group), which 
explicitly demonstrates quantum interference of (coherent) 
states. 



1 INTRODUCTION 



In this paper we consider some starting points in the appli- 
cations of a new numerical-analytical technique which is 
based on the methods of local nonlinear harmonic analy- 
sis (wavelet analysis, generalized coherent states analysis) 
to the quantum/quasiclassical (nonlinear) beam/accelerator 
physics calculations. The reason for this treatment is 
that recently a number of problems appeared in which 
one needs take into account quantum properties of par- 
ticles/beams. We mention only two: diffractive quan- 
tum limits of accelerators (achievable transverse beam spot 
size) and the description of dynamical evolution of high 
density beams by using collective models [1]. Our starting 
point is the general point of view of deformation quantiza- 
tion approach at least on naive MoyalAVeyl/Wigner level 
(from observables to symbols) (part 2). Then we present 
some useful numerical wavelet analysis technique, which 
gives the most sparse representation for two main operators 
(multiplication and differentiating) in any Hilbert space of 
states. Wavelet analysis is a some set of mathematical 
methods, which gives us the possibility to work with well- 
localized bases (Fig.l) in functional spaces and gives for 
the general type of operators (differential, integral, pseu- 
dodifferential) in such bases the maximum sparse forms. 
The approach from this paper is related to our investiga- 
tion of classical nonlinear dynamics of accelerator/beam 
problems [2]-[10]. The common point is that any solu- 
tion which comes from full multiresolution expansion in all 
time scales gives us expansion into a slow part and fast os- 
cillating parts. So, we may move from coarse scales of res- 
olution to the finest one for obtaining more detailed infor- 
mation about our dynamical process. In this way we give 
contribution to our full solution from each scale of resolu- 
tion or each time scale. The same is correct for the contri- 
bution to power spectral density (energy spectrum): we can 
take into account contributions from each level/scale of res- 
olution. Because affine group of translations and dilations 
(or more general group, which acts on the space of solu- 
tions) is inside the approach (in wavelet case), this method 
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Figure 1 : Localized contributions to beam motion. 



2 QUASICLASSICAL EVOLUTION 

Let us consider classical and quantum dynamics in phase 
space £1 = R 2m with coordinates (x, £) and generated by 
Hamiltonian H{x,£) £ C°°(tt;R). If : n — > fl is 
(classical) flow then time evolution of any bounded clas- 
sical observable or symbol b(x,£) 6 C°°(£!,i?) is given 
by h{x,£) = &($^0,£)). Let H = Op w (H) and 
B = Op w (b) are the self-adjoint operators or quantum 
observables in L 2 (R n ), representing the Weyl quantization 
of the symbols TL, b [12] 
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where u G S{R n ) and B t = e iHt/h Be -im/h be the 
Heisenberg observable or quantum evolution of the observ- 
able B under unitary group generated by H. B t solves 
the Heisenberg equation of motion Bt = (i/K) [H, Bt] . Let 
bt(x,£,;Ti) is a symbol of Bt then we have the following 
equation for it 

k = {H,b t } M , (1) 
with the initial condition bo(x,£ i ,h,) = b(x,£ > ). Here 
{/i 9}m(x, £) is the Moyal brackets of the observables 



f,g e C°°{R 2n ), {f,g} M (x,0 = ftg - g$f, where ftg 
is the symbol of the operator product and is presented by 
the composition of the symbols /, g 



(/«<?)(*, 



i<r,p> /h-\-i<uj,T> fh 



■f(x + u>, p + £) • g(x + r, t + ^dpdrdrdu). 

For our problems it is useful that {/, #}m admits the formal 
expansion in powers of h: 

{f,g}M(x,Z)~{f,g} + 2- j ■ 

]T (-l) l/31 • (d?fDPg) ■ (d?gD«f), 

\a+0\=j>l 

where a — (ai, . . . , a n ) is a multi-index, \a\ = at\ + . . . + 
a n , D x = —ihd x . So, evolution (1) for symbol b t {x,^;%) 
is 



h = {H,b t } + ± £ 

|a|+/3|=j>l 

n\d?HD%) ■ {dh t D<* x H). 



(2) 



At ?i = this equation transforms to classical Liouville 
equation 

b t = {H,b t }. (3) 

Equation (2) plays a key role in many quantum (semiclas- 
sical) problem. Our approach to solution of systems (2), 
(3) is based on our technique from [11] and very useful 
linear parametrization for differential operators which we 
present now. Let us consider multiresolution representa- 
tion ... c Vi C Vi C Vq C V-i C V-2 ■ ■ .. Let 
T be an operator T : L 2 (R) -> L 2 (R), with the ker- 
nel K(x,y) and Pj : L 2 (R) Vj (j e Z) is projec- 
tion operators on the subspace Vj corresponding to j level 
of resolution: (Pjf)(x) = J2k < f><Pj,k > <Pj,k{x). 



Let Qj 
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Pj is the projection operator on the 



subspace Wj then we have the following "microscopic or 
telescopic" representation of operator T which takes into 
account contributions from each level of resolution from 
different scales starting with coarsest and ending to finest 
scales [13]: T = V, z {QjTQ 3 + Q J TP J + PjTQj). We 
remember that this is a result of presence of affine group in- 
side this construction. The non-standard form of operator 
representation [13] is a representation of an operator T as 
a chain of triples T = {Aj , Bj , Tj }jez, acting on the sub- 
spaces Vj and Wf Aj : Wj Wj,Bj : Vj Wj,Tj : 
Wj — > Vj, where operators {Aj, Bj,Tj}j e z are defined 
as A 3 = Q 3 TQ 3 , Bj = Q 3 TP 3 , Tj = PjTQj. The 
operator T admits a recursive definition via 



-ij+i 
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where Tj 



PjTP 3 and Tj works on Vj 



Vj ^ Vj. 



It should be noted that operator Aj describes interaction 
on the scale j independently from other scales, opera- 
tors Bj , Tj describe interaction between the scale j and 



all coarser scales, the operator Tj is an "averaged" ver- 
sion of Tj_i. We may compute such non-standard repre- 
sentations of operator d/dx in the wavelet bases by solv- 
ing only the system of linear algebraical equations. Let 
m = J ip(x — £)4-(p(x)dx,£ G Z. Then, the representa- 
tion of d/dx is completely determined by the coefficients 
rg or by representation of d/dx only on the subspace Vq. 
The coefficients m,l e Z satisfy the usual system of lin- 
ear algebraical equations. For the representation of op- 
erator d n /dx n we have the similar reduced linear system 
of equations. Then finally we have for action of operator 
Tj{Tj :Vj — ► Vj) on sufficiently smooth function /: 

(Tjf)(x) = £ (2- J 'X>&*-<) 
where Lpj^ix) = 2~°/ 2 ip(2~ 3 x — k) is wavelet basis and 
fj,k-i = 2- j/2 / f(x)<p(2-*x — k + £)dx 



are wavelet coefficients. So, we have simple linear para- 
metrization of matrix representation of our differential op- 
erator in wavelet basis and of the action of this operator on 
arbitrary vector in our functional space. Then we may use 
such representation in all quasiclassical calculations. 

3 WIGNER TRANSFORM 

According to Weyl transform (observable-symbol) state or 
wave function corresponds to Wigner function, which is 
analog of classical phase-space distribution. If ip(x, t), x 6 
R n satisfies the Schroedinger equation 

-(h 2 /2)AiP + V^ 



ihd t ip 

and W is the Wigner transform of tp 

W{t,x,v)= Je-^jfrx+ihffly)- 
^{t,x-{h/2)y)dy, 



(4) 



(5) 



then W satisfies the pseudo-differential (i/>DO) Wigner 
equation 

d t W + vd x W -{i/K)P(V)W = 0, (6) 
where i/>DO operator P(V) is 

P(V)f(x, u) = j^J e-™y [V(x + *y) - 

V(x-*y)] ■(fe?ytf(x,Od$)dy (7) 

In quasiclassical limit % — > the operator P(V) converges 
to — d x V ■ d v . We consider it in [1 1]. On Fig. 2 we present 
calculations [14] of Wigner transform for beam motion, 
represented by four gaussians, which explicitly demon- 
strates quantum inteference in the phase space. 
We give more details in [11]. 
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Figure 2: Wigner transform 
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